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1. Introduction 

The issue of form invariance symmetries [1] has been invoked recently to extract 
different evolution regimes from Friedmann- Robert son- Walker (FRW) and Bianchi type 
V cosmologies [2]. This is the case when we obtain the so-called phantom cosmologies [3j- 
[6] from ordinary regimes like universes filled with barotropic fluids [7] . These phantom 
models generate solutions that contain negative pressure situations, therefore promoting 
positive accelerated expansions [H],[9]. The case where the gravitational source is a 
bosonic field was investigated in [T],[3],[5], and the question that follows would be what 
happens when the field promoting the universe's expansion is a fermion whose dynamics 
is described by the Dirac equations. This kind of model has been analyzed in several 
works [To], [11] where it was shown that regimes with ordinary matter/dark energy 
transitions are possible solutions of the Einstein and Dirac equations. Taking these 
ideas into account, the main purpose of this work is to investigate the form invariance 
transformations (FIT) that result when the Klein-Gordon dynamics is replaced by the 
Dirac equations in a FRW space-time. The paper is structured as follows: in section 
II we review the basic ideas behind the FIT in cosmology; in section III we focus on 
the scalar field case. In section IV we analyze the symmetries behind the fermionic 
formulations and analyze the cosmological constant group with a fermionic seed fluid. 
In section V we introduce the dual transformation and the associated "phantomization" 
process. Finally in section VI we present our conclusions. 

2. Internal symmetry in FRW 

An interesting method of obtaining new exact solutions of the Einstein equations from 
already existing ones, by using the FIT, has been developed in a series of papers [l]-[6]. 



In spatially fiat perfect fluid FRW cosmologies, such transformations can be viewed 
as a prescription relating the quantities a, H, p and p in a given initial scenario to 
quantities a, p and p corresponding to a new cosmological model. 

As a starting point we present a quick review of this internal symmetry with the 
purpose of investigating it in bosonic and fermionic cases. As is well-known, the Einstein 
equations for a fiat FRW cosmological model with scale factor a, and filled with a perfect 
fluid with energy density p and pressure p, are 



m-m- 



= p, 

p + 3H{p + p) = 0, 




(2) 



where H = a/a. The corresponding FIT is given by [T] 
P = P(P) , 



(3) 




(4) 
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P--p+(f)'^. + p)|. (5) 

where p = p{p) is an arbitrary invertible function and the set of transformations (l3])-(l5l) 
give rise to the form invariance symmetry (FIS) group . We are interested in analyzing 
the perfect fluid case, where heat conduction, shear and bulk viscosities are absent. 
Furthermore, we assume a barotropic equation of state for the pressure which reads 
p = (7 — Hence, it follows that barotropic indices of the original 7 and transformed 
fluid 7 are related by 

For later application it will be useful to investigate the transformation rule for the 
quantities which characterize the geometry and the fluid when the transformation is 
generated by [1] 

p = n^p, (7) 

where n is the constant parameter of the FIS group. Hence, for the geometrical 
quantities we get 

H = nH, a = a^, (8) 

and 

iZ; = 7 - 1 = - - 1, (9) 

n 

p = wp= (J^-l^p, (10) 

for the state parameter w and the pressure of the fluid. 

Let us analyze the equation of state ffTOl) generated by transformation ([7]) when the 
seed fluid represents a normal fluid with 7 > 0. We can infer: (i) when < n < 7 
it follows that w > and the cosmological fluid behaves as a normal fluid; (ii) when 
n > 7 it follows that —l<w<0 and the cosmological fluid behaves as a quintessence 
constituent; (iii) when n < it follows that w < —1 and the cosmological fluid behaves 
as a phantom constituent. 

3. The bosonic case 

In this section we analyze the transformation rules for a scalar fleld 0, self-interacting 
through a potential V{(j)), under the simple transformation generated by (171) [1]. 
Applying Eqs. and ([5]) to the dynamical quantities associated with the scalar fleld, 
i.e., its energy density p^f, = 0^/2 + V{(f>) and pressure = 0^/2 — V{(f)), we obtain 

0' = n0', (11) 

1^(0) =n(n-l)^ + nV(0). (12) 
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Note that t^(0) given in the above equation is not properly a potential, since it depends 
on the kinetic and potential terms of the seed fluid, but it will become a potential 
once the cosmological solution of the seed fluid is known. In other words, the FIT 
acts on the internal space of the Friedmann's cosmology so the cosmological time is not 
affected by those kinds of transformations. The knowledge of the seed solution (j){t) 
allows us to invert it obtaining t = t{(f)), <p{t) = <p{t{(j))) and <p = 4'{4')- Using the 
latter in Eq. ( ITTil we get = n(j) {(j)) and (p = \fn<p + 0o where is an integration 
constant. Then, the kinetic and the potential terms become 0^(0) = — 0o)/v^) 

and V^(0) = 0o)/v^), respectively. Inserting these last two relations into the FIT 

(ITT]) and (fT2!) leads to the transformed potential V^(0) which is expressed as a function 
of the transformed potential 0. 

In order to study the action of the group, we will now generate the set of power 
law solutions starting from a particular seed one. 



3.1. Bosonic power law group 

Let us solve the Einstein-Klein-Gordon equations for a free scalar field, V = The 
solutions of those equations are given by 



aW = (±^)l/^ = ±y-ln|t|, 1/ = 0. (13) 

where the branches a^^"^ and a*-"^ are defined for t > and t < respectively. Here we 
call attention to the fact that the ± branches of the scalar field are not related to the 
branches a^^^ and a^~^ of the cosmic scale factor. Now, we use the latter as a seed set 
for the FIS group. To do that, we first calculate the transformed quantities when the 
group is acting on the above solutions 

02 = n02, (14) 
V{4>)=n{n-l)^. (15) 

Finally combining the seed set f[T^ with the above equations we obtain the power law 
solutions together with the scalar fields and corresponding potentials [5] 

^(±) = (±^)n/3^ (16) 



V{<P) 




;i7) 



For the identity transformation, n = 1, the latter equations (fT6l) -( fT8l) reduce to the 
seed solution (fT3l) . We conclude that power law solutions can be generated from a seed 
solution corresponding to the free scalar field. 
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4. The fermionic case 



In this section we extend the FIT to a fermionic field satisfying the Dirac equation in 
curved space time. As required by equations (1) and (2) we must calculate the energy 
density and the hydrostatic pressure of the fermionic field. We have [H] 

m{i^ip) + V, (19) 



(20) 



where ip and ip = ip'^'-f'^ are the spinor field and its adjoint, respectively. 

To obtain the transformation properties of and we express the potential in 
terms of the scalar invariant X = {ip^Y and the pseudo-scalar invariant Y = {i^p'j^ipY. 
To this end, we assume a generic potential V = V{X,Y), that includes, among others, 
the Nambu-Jona-Lasinio potential [TT]. For this kind of potential the hydrostatic 
pressure becomes 



dV dV 
2X— + 2Y— - V. 



dX dY 
Applying the transformation rules ([3])- 
hydrostatic pressure (12T1) we find that 

/J+\7 = n2(mVX + y), 



(21) 

©, ([7j) to the energy density (fT9|) and the 



-dV -dV 



mV X 



n 



2X 



.dV 
dX 



dV 
2Y— + m 
BY 



(22) 
(23) 



Inserting V, from Eq. fl22|) . into Eq. fl23l) we obtain the following equation 



2^fX dX 



dX dY n 



dV 



dY 
^dX 



-dY 
dY 



Y 



n 



0. 



(24) 



All transformations satisfying the last condition (1241) represent internal symmetries of the 
Einstein equations with a fermionic source, provided that condition becomes an identity 
for any potential and particle mass. In any other case, the above equation becomes a 
partial differential equation whose solutions would lead to particular transformations 
depending on the transformed quantities. We must discard them because a symmetry 
transformation depends on a set of parameters or functions which do not include the 
quantities to be transformed. Hence, a necessary and sufficient condition for Eq. fl2^ to 
be satisfied for any potential and particle mass is that the two terms in square brackets 
must vanish separately, i.e., 

X , , 

(25) 



dX 
dY 

^dX 



dY 
^dY 

dY 



n 
Y 

n 



(26) 
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For a general massive fermionic field whose potential depends only on X, the 
solution of fl25l) along with (l22l) leads to the general FIT 

X = X'^, (27) 
V = n^V + n^mVX - mX'^l'^ , (28) 

where we have normalized to one the multiplicative integration constant in Eq. ( 1271) . 
Because the fermionic masses represent two different particles, we consider both masses 
as free parameters of the model. Then, Eqs. ( 1271) and ( l28i) define a FIS of the Einstein- 
Dirac dynamics. Also, from Eqs. (fT9!) and (I2T!) the equation of conservation ([2]) becomes 

(myX + 2X^ 1 (X + 6Xi/) = 0. (29) 
\ aX ) 

From (29), it is easy to see that the two particular cases {\)V = Vq, m 7^ and (ii) m = 0, 
y 7^ Vo lead to the same solution, X oc a~^, so that V can be treated interchangeably 
as a function of X or a function of a. 



4-1. Fermionic power law group 

Following the same methodology as the bosonic scalar field, let us investigate two simple 
cases, (i) the free massive fermionic particle and (ii) the massless fermionic particle with 
a polynomial potential. 

In the free case, the energy density f|T9l) and pressure fl2T]) of the fermionic particle 
are given by 

p = mVx, p=0, V = 0. (30) 
Solving the Einstein equations we obtain 

„<.).„„(±,V. X^il(^)\ (31) 

where oq is an integration constant. Suppose now that we are interested in comparing 
this cosmological model with another one driven by a potential, which is assumed to 
depend only on X. This can be done by inserting the seed configuration (l30 l) -( l3Ti) into 
the transformations ([5]), ([7]), ([8]) and (l271) - (!28ll . Then, we find the full set of power 
law solutions along with the potential, energy density and pressure characterizing the 
transformed fermionic configuration: 



a" ' = 0-0 



±t:f'''\ a„ = a;, (32) 



V = n^mX^/^" _ mX^/2, (34) 
p = n'mX^/'", p=(--l]p, (35) 



while the Einstein equations for the new configuration are 

r2 



3H' = p, -2H = p + p. (36) 
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The transformed potential V has a vanishing hmit for large cosmological time matching 
smoothly with the seed solution (l30|) - fl3T]) . While the seed solution fl30|) - fl3T]) may be 
identified with some kind of fermionic dust matter the new solution ( I32l) -(l35l) may be 
associated with a fermionic barotropic fluid whose equation of state is 7 = l+p/ p = 1/n. 

The Dirac equations in a FRW background for the seed and transformed 
configurations read 

3 d\'^ 
^^ + -i^7^ + ^m7V + ^7°^ = 0, (37) 
/ dip 

3 d\'^ 
ip + -Hip + zm7V + «7°— = 0, (38) 
2 dip 

along with the corresponding equation for the adjoint ip and ip. In order to preserve 
the characteristics of the fermionic particles, we assume that the Dirac matrices remain 
invariant under a FIS transformation. 

Solving the Dirac equation ( 1381) for the new configuration (I32l) -( l35l) . we obtain the 
spinor field solution 



3/2 



/ bi e"*'^™^" 

^* ^mrriTn 



(39) 



where 



and 



\M^+\h2?-\d,\^-\d2 



( 4 



>-) 

V3m/ 



> 0, 



/ 4 y-" 

V3m/ 



3(l-n)/n 



dt. 



(40) 



(41) 



For n = \ the latter becomes the spinor field corresponding to the seed solution (1301) - 

(EH). 

Now, we consider the massless case and investigate a spinor field driven by a 
potential V = V{X) depending on X. Here, the energy density p and the pressure 
p associated with the spinor field are given by: 



dV 



(42) 



The imposition of the energy density transformation p = n^p reduces the transformation 

rules (EZD-iEHD to 

X = X", V = nW - mX"/2. (43) 
By choosing the seed potential V = VqX", with Vq and a constants, we obtain 

p = p={2a- l)p, (44) 
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and from the Einstein equations it follows 

„(±„V.«, (45, 

In this case power law solutions for the transformed fermionic configuration read 



a = ao(±t)"/'", ao = (46) 




V = nVoX"/" - mX^/\ (48) 

p = nVoX"/", p=^^-{jp. (49) 

In the massless case the solution of the Dirac equation ( l38|) for the latter 
configuration (H^ - fH^ is given by 



where 



and 




In summary, in both the free and massless cases we have found that the Einstein- 
Dirac equations admit a power law group. In the former it is generated by fermionic 
dust matter, while in the latter it is generated by a polynomial potential V cx X". 
In some sense, there is a resemblance between a scalar field driven by an exponential 
potential and a massless Dirac field driven by a polynomial potential. This comparison 
is justified by the fact that both induce equations of state of the form p oc p. Due to 
the importance of the use of power law solutions in cosmological models, it would be 
interesting to investigate the existence of power law groups for other kinds of fields. 



4-2. Cosmological constant group 

In this subsection we consider a mixture of the above two examples namely, the fermionic 
power law group in the free and massless cases summarized in the seeds given by Eqs. 
(l30l) and (H2l) . The former leads to a pressureless dust dominated solution whereas the 
latter leads to a de Sitter solution in the constant potential case. It is well known that 
the ACDM cosmological model interpolates between these two asymptotic behaviors. 
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Thus, using this cosmological model as a seed one 



V - 
X 



A, 
16 
9m2 



p = mV X 

6 



a J 



ao 



±2 



sinh 



3A 



A, 



2/3 



P 



-A, 



/3A 2 

where A is the cosmological constant, the FIS gives 



V = n^A 

p = n^ (mX'/^^ + A' 
p = n{l- n)mXi/2n _ ^2^^ 
X ' - ^ ' 



X, 



V9m2y 



ao 



±2 



3A 



sinh 



3A 



2n/3 



(53) 
(54) 

(55) 

(56) 
(57) 
(58) 

(59) 
(60) 



Asymptotically, in the limit t oo, the new potential tends toV—* n^A for n > 0. 
Therefore, the last set of equations shows that there exists a large set of models which 
generalize the ACDM one, i.e., the class of model generated by the set of FIT ([7]) with 
cosmological constant A ^ n^A. They describe universes containing a perfect fluid and 
with cosmological constant evolving as a power-law dominated phase at early times, 
where the scale factor behaves as a oc t^"/^. Finally, such universes end in a stable 
de Sitter accelerated expansion scenario with a scale factor behaving as a oc e"v^^*, 
representing a possible dark energy final stage. 

Finally, for the new configuration (!32|) - (l35|) the spinor field solution of the Dirac's 
equation (l38i) becomes 



where 



and 



3/2 



^* ^mmTn 



(61) 



W+\h2?-\d,\^-\d2 



{-X 

\3mJ 



> 0, 



/ 4 y-n 
\3mJ 



3(l-n)/n 



dt. 



(62) 



(63) 



This spinor depends on the cosmological constant through the transformed scale factor 
d. 
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5. Duality and Phantom cosmologies 



The "phantomization" process means that the energy density of the expanding 
cosmological model being analysed must satisfy p > 0, or equivalently the weak energy 
condition (WEC) must be violated, so that (p+p) < [3]. In terms of the Hubble factor 
H, these two conditions become H > and H > 0. Therefore, one can infer two cases, 
namely, (a) p has an asymptote p — > A for t — oo, or (b) p grows without limit. In the 
first case the scale factor tends to a de Sitter solution a — > exp ^ A/3t. In the second 
case if it is assumed that the energy density has the asymptotic behavior p p^a^ with 
/c > 0, the asymptotic solution of the Friedmann equation reads 

2v/3 



,(-) 



,(+) 




-| 2/k 



t < u 



t > u 



(64) 



(65) 



k./p^{t - toj 

The expanding solution a^~^ is defined for t < to and ends in a big rip at t = to? since the 
scale factor diverges at the finite time to resulting in a future singularity. By contrast, 
the contracting solution a^~^^ begins from a past singularity at t = to- To sum up, a^~^ 
is the "phantomization" of the solution l/a^~^ which ends in a big crunch at t = to- In 
terms of FIS the "phantomization" process means that there is a duality between two 
solutions of the Friedmann equation, i.e., between l/a*-"-* and a*-"-* [12]-[16]. In fact, the 
class of FIT generated by n = — 1: 

P^P = P, P + P ^ P + P = -{? + p), (66) 

with 



H = -H, 



1/a, 



(67) 



transforms a contracting scale factor, H < 0, satisfying the WEC, p + p > into 
an expanding one, H > 0, which violates the WEC, so that p + p < 0. The latter, 
dubbed the dual transformation, is crucial because of its applicability as a method of 
transforming a conventional cosmological model into a phantom one by performing a 
FIT. 

In the bosonic case the expressions (ITT]) and (IT^ become 

0' = -0', (68) 
V^(0) =02 + v^(0). (69) 

and this implies that (p = icj). That is, the transformed scalar field is related to the 
original one by a Wick rotation, and we finally have 



P 



(70) 



The sign in the kinetic part of the energy density, 0^ = —0^ < 0, indicates that what 
we actually have now is a phantom cosmology, with a phantom field 0, driven by the 
potential t^(0), which is a real function of 0. 
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For the fermionic matter, in the case where the potential depends on X, the 
"phantomization" induces the following transformation rules for the variables X and 
V 

X ^ i (71) 
V = mVx (l - ^) + V- (72) 

and 

7 = -7, (73) 
for the barotropic index of the spinor field. Hence, there is a duality between a 
contracting universe filled with an ordinary spinor field with 7 > 0, and an expanding 
universe filled with a phantom spinor field with 7 < 0. 



6. Conclusion 



We have found the symmetry transformations under which the Einstein equations, for 
a spatially fiat FRW space-time filled with bosonic or spinor fields, preserve their form. 
This group of transformations has been used to obtain power-law solutions from a seed 
one, for bosonic or fermionic fluids. We have given special attention to the n = —1 
case, which expresses the duality between non-accelerated and accelerated scenarios 
and vice versa. Thus, starting from a contracting spatially fiat FRW cosmological 
model we get, after using the dual transformation, a super-accelerated spatially fiat 
FRW cosmological model, i.e., the "phantomization" of the model. We have shown 
that bosonic and fermionic fields behave differently under a dual transformation, the 
former becomes imaginary, whereas the latter changes the sign of its phase and inverts 
its asymptotic limits. Here we have extended the analysis of the bosonic case, previously 
studied in the works [3] and [5] , by discussing the role of the seed fluid, the duality and 
"phantomization" of the solutions. The analysis of the cosmological constant group 
with a fermionic seed fluid indicates that one can obtain a model for a universe which 
evolves from a power-law dominated phase to a stable de Sitter accelerated expansion. 
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